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Hidden connections between regression models of wind tunnel strain gage 
balance calibration data are investigated. These connections become visible 
whenever balance calibration data is supplied in its design format and both 
the Iterative and Non-Iterative Method are used to process the data. First, 
it is shown how the regression coefficients of the fitted balance loads of a force 
balance can be approximated by using the corresponding regression coefficients 
of the fitted strain gage outputs. Then, data from the manual calibration of the 
Ames MK40 six— component force balance is chosen to illustrate how estimates 
of the regression coefficients of the fitted balance loads can be obtained from the 
regression coefficients of the fitted strain gage outputs. The study illustrates 
that load predictions obtained by applying the Iterative or the Non Iterative 
Method originate from two related regression solutions of the balance calibration 
data as long as balance loads are given in the design format of the balance, gage 
outputs behave highly linear, strict statistical quality metrics are used to assess 
regression models of the data, and regression model term combinations of the 
fitted loads and gage outputs can be obtained by a simple variable exchange. 
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= regression coefficients of the forward normal force gage output ( Iterative Method ) 
= axial force 

= regression coefficients of the aft normal force gage output ( Iterative Method) 

= regression coefficients of the forward side force gage output ( Iterative Method ) 

= square matrix; used by the load iteration process 

= rectangular matrix; used by the load iteration process 

= regression coefficients of the aft side force gage output ( Iterative Method) 

= regression coefficients of the rolling moment gage output ( Iterative Method) 

= regression coefficients of the axial force gage output (It erative Method) 

= part of matrix G that contains loads 
= load matrix 

= part of matrix G that contains absolute value and non-linear terms 
= load iteration step index 

= normal force at the forward normal force gage of the balance 
= normal force at the aft normal force gage of the balance 
= electrical outputs of the forward normal force gage 
= electrical outputs of the aft normal force gage 
= electrical outputs of the forward side force gage 
= electrical outputs of the aft side force gage 
= electrical outputs of the rolling moment gage 
= electrical outputs of the axial force gage 
= rolling moment 
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= side force at the forward side force gage of the balance 
= side force at the aft side force gage of the balance 

= regression coefficients of the aft normal force ( Non-Iterative Method) 

= regression coefficients of the forward side force ( Non-Iterative Method) 

= regression coefficients of the aft side force ( Non-Iterative Method) 

= delta gage output vector or matrix 

= regression coefficients of the axial force ( Non-Iterative Method) 

= regression coefficients of the forward normal force ( Non-Iterative Method) 
= regression coefficients of the rolling moment ( Non-Iterative Method) 

= index of regression coefficient of a dominant regression model term 
= index of regression coefficient of a dominant regression model term 
= index of regression coefficient of a dominant regression model term 
= index of regression coefficient of a dominant regression model term 
= index of regression coefficient of a dominant regression model term 
= index of regression coefficient of a dominant regression model term 

I. Introduction 


Different analysis approaches are used in the wind tunnel testing community to predict balance loads 
from measured strain-gage outputs during a wind tunnel test. One group of analysts, for example, processes 
balance calibration data by first fitting strain-gage outputs of the balance as a function of the applied balance 
loads. In that case, an iteration scheme is needed so that balance loads can be predicted from measured 
strain-gage outputs during a wind tunnel test. This analysis approach is called the Iterative Method (see 
Refs. [1], [2], and [3] for more detail). 

In principle, the balance calibration experiment defines loads as independent variables and gage outputs 
as dependent variables. Some analysts prefer to switch the independent and dependent variables that the 
balance calibration experiment defines. This alternate analysis approach is called the Non Iterative Method. 
In this case, no iteration is needed to predict loads from gage outputs during a wind tunnel test because the 
applied loads are directly fitted as a function of measured gage outputs (see Ref. [4] for more details) . 

Regression models used by the Iterative Method and Non-Iterative Method are derived from the same 
balance calibration data set. Therefore, in theory, they should contain the same information about the 
behavior of the balance even though, in one case, gage outputs are fitted as a function of balance loads 
and, in another case, balance loads are fitted as a function of gage outputs. The present paper studies the 
relationship between the regression models of the two fundamentally different analysis approaches in more 
detail as the balance characteristics themself must be contained in (i) math terms that are selected for the 
regression analysis and (ii) the sign and magnitude of the regression coefficients. 

In general, it is recommended to process a balance calibration data set in its “design” load format. In 
other words - a force balance should be analyzed in force balance format, or, a moment balance should be 
analyzed in moment balance format. Then, the primary sensitivities of all gages of the balance exist (see 
Ref. [5] for more details). This characteristic also means that, in an ideal case, the gage outputs of single 
gage loadings are located along a straight line when plotted versus the corresponding single gage loads (a 
single gage loading is a setup during the calibration of a balance that applies a single load component to the 
balance while simultaneously keeping the magnitude of all other load components close to zero). The gage 
outputs of the remaining combined loadings will also be in the vicinity of this straight line. 

Figure 1, for example, shows the output of the forward normal force gage of a force balance plotted 
versus the forward normal force. Both single gage loadings and combined loadings are depicted. It can be 
seen that the gage output is more or less proportional to the corresponding primary gage load. The required 
constant of proportionality is the inverse of the primary gage sensitivity of the gage. Therefore, a math 
model term combination selected to fit gage outputs as a function of balance loads could also be used to 
approximate the fit the corresponding primary gage load as a function of the gage outputs (and vice versa). 
It is only required to switch primary loads and gage outputs in the related regression models. 

Now, the following question emerges: Can a direct connection between the coefficients used by the 
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regression model of the Iterative and Non-Iterative Method be established if the primary loads and gage 
outputs are switched in the related regression models? Let us assume that the answer to the question is 
“yes.” This means that the final balance load prediction accuracies of the Iterative and the Non-Iterative 
Method are more closely linked than a casual observer would suspect. It would also help users of the Iterative 
Method to gain confidence in the load prediction accuracy of the Non-Iterative Method (and vice versa). 

First, in order to find an answer to the question posed in the previous paragraph, basic elements of 
the Iterative and Non-Iterative Method are reviewed for a typical force balance. Then, using a calibration 
data set of a force balance as an example, the connection between the regression coefficient sets will be 
investigated in more detail. 


II. Balance Calibration Data Analysis Methods 

A. Iterative Method 

This section describes the balance calibration data analysis assuming the Iterative Method is applied 
to a force balance. Basic elements of the method are discussed in great detail in Ref. [1]. Therefore, only 
an abbreviated description of the application of the method to a force balance is given in this section. In 
principle, the Iterative Method is a two step process. First, gage outputs are fitted as a function of calibration 
loads. Then, the regression coefficients of the gage outputs are used to construct a load iteration process so 
that balance loads can be predicted from measured gage outputs during a wind tunnel test. 

Data from the calibration of a force balance may be used to illustrate the application of the Iterative 
Method. It is assumed that (i) data of a six-component force balance is analyzed and that (ii) the loads are 
given in force balance format. Therefore, the regression models of the six gage outputs can be expressed as 
a function of the balance loads using the following equations: 
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The above equations highlight the fact that the regression model of each gage will be dominated by the 
influence of the primary gage load. Now, the six balance loads need to be computed iteratively after the 
completion of the regression analysis. The following iteration equation in combination with a load iteration 
process may be used for that purpose (see Ref. [1] for a description of the iteration process): 


F i = 

c^ar 


C 1 1 C 2 

H ;i 
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changes for each iteration step 


(2) 


Equation (2) is a matrix equation. Two matrices used in Eq. (2), i.e., Ci and C 2 , are derived from 
the regression coefficients of the gage outputs that are defined in Eqs. (la) to (1/). The vector AR has 
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the gage outputs that are measured when the balance experiences a load. Matrix H is constructed from the 
intermediate load estimates of the previous iteration step. The load iterations typically converge after 5 to 
10 iterations assuming that a tolerance of 0.0001 % of capacity is used to test for convergence. 

B. Non-Iterative Method 

The balance calibration data of a force balance may also be analyzed using the Non-Iterative Method. 
Differences between the Non-Iterative Method and the Iterative Method are discussed in great detail in 
Ref. [4]. Therefore, only basic elements of the Non-Iterative Method are reviewed in this section. 

In principle, the Non-Iterative Method exchanges the independent and dependent variables that the 
Iterative Method uses. Now, gage outputs become independent variables and balance loads become dependent 
variables as far as the regression analysis of the balance calibration data is concerned. The Non-Iterative 
Method is a one step process as the loads are directly fitted as a function of the measured gage outputs. 
Consequently, no load iteration is required to predict loads from gage outputs during a wind tunnel test. 

Again, data of a six-component force balance may be used to illustrate the application of the Non- 
Iterative Method. It is assumed that the calibration data of the force balance is given in force balance format. 
Then, we get the following six regression models for the analysis of the balance calibration data: 


N 1 
N2 

51 

52 
RM 
AF 


J7o^ 

intercept 

intercept 

intercept 

intercept 

intercept 


Co 


intercept 


+ r/£ ■ R1 + 

dominant 

+ f 3p • R2 + 

dominant 

+ '7cr • R3 + 

dominant 

+ 8(f) • RA + 

dominant 

+ fl^fj * Rb + 

dominant 

+ Co> * ^6 + 

dominant 


(за) 

( зб ) 

(3c) 

(3d) 

(3e) 

( 3 /) 


The above equations highlight the fact that the regression model of each load component will be 
dominated by the primary gage output. Loads are computed during a wind tunnel test by using the measured 
gage outputs as input for the regression models of the loads that are defined in Eqs. (3a) to (3/). 

In general, the Non-Iterative Method has the advantage that it is a one-step method. No iteration 
is needed to compute loads from measured strain-gage outputs during a wind tunnel test. An analyst, 
however, must not forget that the Non-Iterative Method ignores the fact that the balance loads are the 
“true” independent variables of the calibration experiment as loads are “applied” and strain-gage outputs 
are “measured” during the calibration of a balance. Therefore, the success of the Non-Iterative Method hinges 
on the fundamental assumption that a switch of the independent and dependent variables of the calibration 
data set does not negatively influence the mathematical description of the “true” physical behavior of the 
balance. In addition, the robustness and reliability of the regression model of each balance load depends on 
the fact that (i) the model does not have near-linear dependencies between terms and that (ii) it consists 
of statistically significant terms (see Ref. [6] for a discussion of these issues). These two requirements also 
apply to regression models of the gage outputs that the Iterative Method uses. 

In the next section of the paper data from the calibration of NASA’s MK40 balance is used to illustrate 
the connection between the regression coefficients of the Iterative Method and Non-Iterative Method. 
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III. Hidden Connection between Regression Models 


Data from the calibration of the NASA Ames MK40 force balance was chosen to illustrate the “hidden” 
connections that exist between the regression coefficients of the strain-gage outputs and the regression 
coefficients of the balance loads. The Ames MK40 balance was manufactured by the Task Corporation. It is 
a six-component force balance that measures five forces and one moment (Al, A 2, 51, 52, AF, RM). The 
balance has a diameter of 2.5 inches and a total length of 17.31 inches. Table 1 shows the load capacity of 
each load component. 

Table 1: Load capacities of the NASA Ames 2.5in MK40 balance. 



Al, lbs 

A2. lbs 

51, lbs 

52, lbs 

RM, in-lbs 

AF, lbs 

CAPACITY 

3500 

3500 

2500 

2500 

8000 

400 


The original balance calibration was performed as a manual calibration. A total of 164 data points 
were taken in 16 load series. The analysis of the balance calibration data was done in several steps. First, 
the given calibration loads were tare corrected for the weight of the balance shell, calibration body, and 
other calibration fixtures. Then, the calibration data was analyzed using both the Iterative Method and the 
Non-Iterative Method that were described in previous sections. 

For simplicity, it was decided to focus on hidden connections between regression coefficients of the 
forward normal force component and regression coefficients of the forward normal force gage output. Con- 
nections between regression coefficients of other load components and gage outputs can be investigated in a 
similar manner. 

First, the Non-Iterative Method was applied to the calibration data of the MK40 balance. The analysis 
started by applying a regression model optimization process to the data so that the regression model of the 
forward normal force component would satisfy a set of widely accepted statistical quality requirements (see 
Refs. [7] and [ 8 ] for a description of the optimization process and the quality requirements). The optimization 
process chose the following 13-term regression model for the forward normal force component: 

N1 = ? 7 o + r/i ■ R1 + ?? 2 • 52 + 773 • 53 + 774 • 55 + 775 - |51| 

+ 776 • |52| + 777 • |53| + rjs ■ |55| + 779 • 55 • 55 (4) 

+ 7710 • 51 • |51| + 7711 • 52 • |52| + 7712 • 55 • |55| 

The symbols 770 , 771 , ..., 7712 are the coefficients of the regression model of the forward normal force 
component. Figure 2a shows the Analysis of Variance result for the regression model of 7VT. The math 
term and f-statistic columns are highlighted using blue rectangles. The f-statistic results confirm that the 
forward normal force gage output 51 is the dominating term in the regression model of the forward normal 
force component A1 as the term 51 has the largest f-statistic magnitude (+982). Figure 2b shows the 

regression model and all coefficients of the forward normal force component. This is the regression model 

that the Non-Iterative Method ultimately uses for the prediction of the forward normal force component 
from the measured strain-gage outputs. 

Now, the Iterative Method was applied to the calibration data of the MK40 balance. Similarly, the 
regression model optimization process identified the following 13-term regression model for the forward 
normal force gage output: 

51 — oq + ai • A1 + a 2 * A2 + <73 • 51 + 04 . * RAI + a§ • |A1| 

+ a 6 • |A2| + o 7 -|51| + a 8 • \RM\ + a 9 ■ RM ■ RM (5) 

+ a i0 • A1 • |A1| + an • A2 • |A2| + a 12 ■ RM ■ \RM\ 

It is interesting to point out that the optimized regression model term combination listed on the right 
hand side of Eq. (5) can also be obtained by simply exchanging each primary gage output with the corre- 
sponding primary gage load in Eq. (4). The symbols a 9 , ai, ..., ai 2 are the regression coefficients of the 
forward normal force gage output 51. Figure 3a shows the Analysis of Variance result for the regression 
model of 51. The math term and f-statistic columns are highlighted using blue rectangles. The f-statistic 
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results confirm that the forward normal force component TV 1 is the dominating term in the regression model 
of AI as the term TV 1 has the largest ^-statistic magnitude (+891). Figure 3b shows the optimized regres- 
sion models and all coefficients of the six gage outputs of the MK40 balance. These coefficients are used to 
assemble the matrices Ci and C 2 that the iteration equation and iteration process of the Iterative Method 
need for the prediction of balance loads from measured strain-gage outputs. 

Now, a question comes up: Can the coefficients defined by the regression model of the forward normal 
force gage output R1 (see Eq. (5)) be used to “reverse engineer” the coefficients of the regression model of 
the forward normal force component TV 1 (see Eq. (4))? An answer to this question can be found by first 
solving the regression model of the forward normal force gage output Rl, i.e., Eq. (5), for the forward normal 
force component TV1. Then, after some algebra, we get the following equation from Eq. (5): 

TV1 = — — + — Rl + — — -TV2 + ^--Sl + — — • RM + — — ■ |TV1| 
a\ a\ a\ a\ a\ a\ 

+ -^--|TV2| + -^--|51| + — — . \rm\ + AM -AM (6) 

ai a 1 a\ a\ 

+ ~ ai ° • TVl • | ATI | + ~ an • TV2 • |TV2| + ~ Ql2 • RM ■ \RM\ 

CL 1 CL\ CL\ 

Unfortunately, Eq. (6) still has the loads TVl, TV2, 51, and RM on the right hand side of the equation. 
These loads, however, can be approximated by using the observations that (i) the original regression models 
of the strain-gage outputs each have a single dominant term and that (ii) all other terms of the regression 
model of the gage outputs are small when compared with the corresponding dominant term. Therefore, after 


simplifying Eqs. (la), (16), (lc), and (le), and inspecting the coefficients given in Fig. 3b, we get: 

AI « ■ TVl where f = 1 

(7a) 

R2 s 

a b p ■ TV2 
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P = 

2 
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A3 ? 

a c a • 51 

where 

a = 
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A5 s 
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The approximations given above can be solved for the four remaining load components that still need 
to be substituted on the right hand side of Eq. (6). Then, we get: 


TVl « 

[ 1 / ai ] • Rl 

where 

[ 1 / ai ] 

= primary seirsitivity of Rl 

(8a) 

TV2 « 

[ 1 / 6 2 ] • R2 

where 

[ 1 / b 2 } 

= primary sensitivity of R2 

(86) 

51 « 

[ 1 / c 3 ] • A3 

where 

[ 1 / C3 ] 

= primary sensitivity of A3 

(8c) 

RM « [ 1 / 64 ] • A5 where [ 1 / e± ] = primary sensitivity of R5 (8d) 

Finally, after using Eqs. (8a) to (8 d) to subsitute the remaining load components on the right hand side 


of Eq. (6) and after some algebra, we get the following approximation of the regression model of the forward 
normal force component TVl of the MK40 balance: 
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The exact solutions of the regression coefficients of the forward normal force component N 1 are defined 
in Eq. (4) and listed in Fig. 2b. These coefficients, i.e. , ... 

Vo , Vi , V 2 , V3 , V4 , ' ' ' 

may be compared with the “reverse engineered’' approximations that were obtained from the regression 
model of gage output Rl. They are defined in Eq. (9) as follows: 

— a 0 1 — a 2 -a 3 - 0,4 

cq ’ oq cq 62 ’ ai C3 ’ cq e4 ’ 

The regression coefficients of the fitted gage outputs are listed in Fig. 3b. Therefore, it is possible to 
compute the “reverse engineered” approximations of the regression coefficients of the forward normal force 
and compare them with the exact solutions that are depicted in Fig. 2b. 

The table in Fig. 4 shows the result of the comparison of the exact and approximated coefficient sets. 
Two observations can be made after inspecting the table: (i) signs of the exact and approximated coefficients 
match; (ii) magnitudes of the exact and approximated coefficients show good agreement. 

The influence of each individual coefficient on the regression model was also investigated in more detail. 
Therefore, the exact solution was modified by replacing one coefficient at a time by its approximation. Then, 
the standard deviation of the load residuals was computed for the modified regression model that consisted 
of one approximated and twelve exact coefficients. The computed standard deviations are shown in the last 
column of Fig. 4. We observe that the largest standard deviation is reported for the case when the coefficient 
of the most significant term, indicated by the t-statistic of +982, is replaced by its approximation. This 
observation is expected as the standard deviation of the load residuals must reach a maximum when the 
most significant coefficient is replaced by its approximation. 

Finally, it is of interest to compare the standard deviation of the residuals of the forward normal force 
component for the three different analysis options that are discussed in the present paper. Table 2 below 
lists the standard deviation as a percentage of the capacity of the forward normal force component. 

Table 2: Standard deviation of the load residuals of the forward normal force N 1. 


BALANCE CALIBRATION DATA ANALYSIS APPROACH 

STANDARD DEVIATION 

ITERATIVE METHOD +=+ exact solution, Eq. (2) 

0.0515 % CAP. 

NON-ITERATIVE METHOD +=+ exact solution, Eq. (4) 

0.0489 % CAP. 

NON-ITERATIVE METHOD ++> approximation, Eq. (9) 

0.0646 % CAP. 


As expected, the standard deviations of the exact solution for the Iterative Method and the exact 
solution of the Non-Iterative Method show excellent agreement as the difference between the two standard 
deviations is only 0.0026 %. The standard deviation of the approximation of the Non-Iterative Method shows 
reasonable agreement with the corresponding exact solution as the difference is 0.0157 %. At this point it is 
important to emphasize that the approximation of the Non-Iterative Method should never be used instead of 
the corresponding exact solution for the calculation of balance loads. The approximation was only developed 
for the present study to show that a “hidden” connection between the regression coefficients of the Iterative 
and Non-Iterative Method can be established. 

IV. Conclusions 

The present study illustrates that “hidden” connections between regression coefficient sets used by 
the Iterative and Non-Iterative Method may exist. It is possible to estimate sign and magnitude of the 
regression coefficients of the fitted loads by using the regression coefficients of the fitted gage outputs as 
long as (i) balance data is analyzed in its design format, (ii) the regression models of the loads and gage 
outputs meet rigorous statistical quality requirements, and (iii) the regression model terms can be obtained by 
simply switching primary loads and gage outputs. Numerical differences between the exact and approximated 
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regression coefficient sets remain. They reflect the fact that the regression models used by the Iterative and 
Non-Iterative Method are still independent regression solutions of a given balance calibration data set. 

Results of the present study may help users of both the Iterative and Non-Iterative Method to better 
understand how the two balance calibration data analysis approaches are connected to each other. It must 
also not be forgotten that the application of the Non-Iterative Method requires a switch of the independent 
and dependent variables that the balance calibration experiment defines. Questions about the validity of this 
variable exchange may come up at some point in time. Then, the approach used in the present investigations 
could be used to assess the influence of this variable exchange on the balance load estimates that are obtained 
by applying the Non-Iterative Method. 
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Fig. 1 Relationship between primary gage output R1 and tare corrected primary gage load N 1. 
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MATH TERM CONTRIBUTES SIGNIFICANTLY IF THE T-STATISTIC OF ITS COEFFICIENT IS GREATER THAN 1.976 

(SIGNIFICANCE = 2.50 %, DEGREES OF FREEDOM OF RESIDUAL =151) 
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TERM 

TERM 

COEFFICIENT 

STANDARD 

T-STATISTIC OF 

P -VALUE OF 

VI F 

VI F 

INDEX 

NAME 

VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 

INTERCEPT 

+ 56.5015 

+0.2313 

+ 244.3080 

- 

- 

- 

2 

R1 

+ 2.9148 

+0.0030 

+981.7952 

< 0.0001 

+2.1407 

[ + 13.3766] 

3 

R2 

+0.0800 

+0.0027 

+29.4708 

< 0.0001 

+2.3276 

[ + 12.5974] 

4 

R3 

+0.0024 

+0.0006 

+4.1644 

< 0.0001 

+ 1.0040 

+ 1.0040 

6 

R5 

-0.0174 

+0.0015 

-1 1.3326 

< 0.0001 

+ 1.5522 

+ 8.7920 

8 

IR 1 1 

-0.0256 

+0.0010 

-26.6581 

< 0.0001 

+ 1.0461 

+ 1.0615 

9 

IR2I 

-0.0157 

+0.0009 

-17.8669 

< 0.0001 

+ 1.0502 

+ 1.0509 

10 

IR3I 

+0.0083 

+0.0007 

+ 12.6431 

< 0.0001 

+ 1.0969 

+ 1.0969 

12 

IR5I 

+0.0250 

+0.0018 

+ 14.0078 

< 0.0001 

+9.8290 

+9.8242 

18 

R5*R5 

— 8.9250e-06 

+ 1.7078e-06 

-5.2261 

< 0.0001 

+9.3529 

+9.3686 

20 

R1*IR1 

— 3.2931 e-05 

+5.0333e-06 

-6.5425 

< 0.0001 

+2.4395 

[ + 13.1262] 

21 

R2*IR2 

—3. 1 870e-05 

+4.1 1 27e-06 

-7.7493 

< 0.0001 

+2.2598 

[ + 1 1.6368] 

24 

R5*IR5 

+8.6645e-06 

+ 1.5918e-06 

+ 5.4431 

< 0.0001 

+ 1.5437 

+8.8223 


Fig. 2a Non-Iterative Method: Analysis of Variance results for regression model of forward normal force N 1. 
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INDEX 

TERM 


N 1 

1 

INTERCEPT 

fio= 

+ 5.6501 50e+01 

2 

Rl 

%= 

+ 2.9 1 4783e+00 

3 

R2 

ti 2 = 

+ 7.997201e-02 

4 

R3 


+ 2.409364e-03 

6 

R5 


- 1 .741 967e-02 

8 

IR1 1 

Tls= 

-2.560737e-02 

9 

IR2I 

fi6 = 

— 1 .566803e-02 

10 

IR3I 

11 7 = 

+8.287214e-03 

12 

IR5I 

n 8 = 

+ 2.498232e-02 

18 

R5*R5 

fi 9 = 

-8.924993e-06 

20 

R1.IR1 1 

fiio= 

-3.293052e-05 

21 

R2*IR2I 

fill= 

-3.1 87047e-05 

24 

R5»lR5l 

fil2 = 

+8.664467e-06 


Fig. 2b Non-Iterative Method: Coefficients of optimized regression model of forward normal force IV 1. 



SOURCE OF 

SUM OF PRESS DEGREES OF MEAN 

F-VALUE OF 

P -VALUE OF 



VARIATION 

SQUARES STATISTIC FREEDOM SQUARE REGRESSION 

REGRESSION 



REGRESSION 

4.8039e+06 

12 

400326.6413 991674.6884 

< 0.0001 



RESIDUAL 

60.9568 81 

.4320 151 

0.4037 

- 



TOTAL 

4.8040e+06 

163 

- 

- 





R-SQUARE 

ADJ. R-SQUARE 

PRESS R-SQUARE 





0.999987 

0.999986 

0.999983 




MATH TERM CONTRIBUTES SIGNIFICANTLY IF THE T-STATISTIC OF ITS COEFFICIENT IS GREATER THAN 1.976 





(SIGNIFICANCE = 

2.50 %, DEGREES OF FREEDOM OF RESIDUAL =151) 




DIAGNOSTIC DESCRIBING NEAR-LINEAR DEPENDENCIES (MULTICOLLINEARITY) BETWEEN MATH 

MODEL TERMS 




MAXIMUM OF VARIANCE INFLATION FACTOR (PRIMARY VALUE) = 9.18 < SELECTED LIMIT OF 10 


PHYSICAL VARIABLES & UNITS - REGRESSION COEFFICIENT ESTIMATES AND STATISTICAL METRICS (Rl) 





REGRESSION MODEL HIERARCHY CHARACTERISTICS = HIERARCHICAL 



TERM 


TERM 


COEFFICIENT 

STANDARD 

T-STATISTIC OF 

P -VALUE OF 

VI F 

VI F 

INDEX 


NAME 


VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 


INTERCEPT 

-19.4655 

+0.0733 

-265.4232 

- 

- 

- 

2 


N 1 


+ 0.3438 

+0.0004 

+890.7522 

< 0.0001 

+ 2.8244 

[+14.9350] 

3 


N2 


-0.0103 

+0.0004 

-26.6471 

< 0.0001 

+ 2.8340 

[+14.9156] 

4 


SI 


-0.0005 

+0.0001 

-5.1875 

< 0.0001 

+ 1.0000 

+ 1.0000 

6 


RM 


+ 0.0010 

+0.0001 

+ 1 1.5307 

< 0.0001 

+ 1.6149 

+ 8.7907 

8 


INI 


+ 0.0029 

+0.0001 

+ 25.4357 

< 0.0001 

+ 1.0467 

+ 1.0467 

9 


IN2I 


+0.0018 

+0.0001 

+ 15.6930 

< 0.0001 

+ 1.0458 

+ 1.0458 

10 


IS1 


-0.0013 

+0.0001 

-1 1.6389 

< 0.0001 

+ 1.1207 

+ 1.1207 

12 


RM 


-0.0013 

+0.0001 

-13.3925 

< 0.0001 

+9.1771 

+9.1771 

18 


RM*RM 


+ 7.529 1 e-08 

+ 1 .57 1 2e— 08 

+4.7919 

< 0.0001 

+8.6760 

+8.6760 

20 


n i * In i 


+ 1.5222e-06 

+ 2.1 583e— 07 

+ 7.0526 

< 0.0001 

+ 2.51 12 

[+13.1553] 

21 


N2*IN2 


+ 1 .5540e-06 

+ 2.1 704e-07 

+ 7.1598 

< 0.0001 

+ 2.4990 

[+13.1 181] 

24 


RM«|RM 


-8.6090e-08 

+ 1.51 96e— 08 

-5.6654 

< 0.0001 

+ 1.6149 

+8.7907 












Fig. 3a Iterative Method: Analysis of Variance results for regression model of gage output Rl. 
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INDEX 

TERM 


Rl 

R2 U 

R3 

R4 

R5 

R6 

1 

INTERCEPT 

a o= 

-1.946546e+01 

-2.1 16026e+o/ 

2 C 

-1.571845e+o/ 

3 

-5.743596e+01 

p 

-1.976257e+0 / 

-1.384679e+02 

2 

n i a 

1=> 

+3.437529e-01 

-4.302222e-\fe 

+6.338849e-/3 

0 

0 / 

+4.835826e-03 

3 

N2 

a 2 = 

-1.032882e-02 

+3.725837e-01 

+ 1.123877e>B3 

+3.8579 18e-03 

-4.199624e-|4 

-2.833567e-03 

4 

SI 

a 3= 

-5.235353e-04 

+ 2.862075e-03 

+4.905652e-01 

-2.1 70326e-02 

+4.852466e-|D4 

-2.901 236e-03 

5 

S2 


0 

+2.409149e-03 

-6.977922e-03 

+5.014106e-01 

-3.528204^t>4 

+5.446800e-03 

6 

RM 

a 4 = 

+ 1.035526e-03 

-4.833653e-04 

+4.792302e-03 

+ 5.291 577e-03 

+ 1.638645e-01 

+2.991 91 2e-03 

7 

AF 


0 

-2.708377e-03 

0 

0 

+ 1.3101 20e-02 

+3.569037e+00 

8 

INI 1 

a 5 = 

+ 2.888430e-03 

+ 1 .4581 lOe— 03 

-2.1 73966e-03 

0 

+3.078684e-04 

+ 1.276945e-03 

9 

IN2I 

a 6 = 

+ 1.786223e-03 

+4.314710e-03 

0 

- 1 .371 239e-03 

0 

+ 1 .4401 65e— 03 

10 

IS1 1 

a 7 = 

-1. 34021 6e-03 

-6.562695e-04 

+ 1 .1 05774e— 02 

+ 1 ,303598e-03 

0 

0 

1 1 

IS2I 


0 

—7.61 1 073e— 04 

+ 1 ,495335e— 03 

+ 1.321 102e-02 

0 

0 

12 

IRMI 

a 8 = 

-1.31 6052e-03 

-1.11 5338e-03 

+3.61 7666e-03 

+ 1.803646e-03 

+ 1 .1 55335e-04 

0 

13 

lAFl 


0 

0 

0 

0 

0 

0 

14 

N1*N1 


0 

0 

0 

0 

0 

0 

15 

N2*N2 


0 

0 

0 

0 

0 

0 

16 

S1*S1 


0 

0 

0 

0 

0 

0 

17 

S2*S2 


0 

0 

0 

0 

0 

-9.028354e-07 

18 

RM*RM 

II 

CTl 

+ 7.529071e-08 

+5.54681 2e-08 

-3.001 823e-07 

- 1.347943e-07 

0 

+2.522531e-07 

19 

AF*AF 


0 

0 

0 

0 

0 

0 

20 

n i * In 1 1 

a io 

+ 1 .5221 83e-06 

+ 1.368479e-06 

0 

0 

0 

0 

21 

N2*lN2l 

a ll= 

+ 1.553984e-06 

+ 1.580928e-06 

0 

-1.270365e-06 

0 

0 

22 

S1*lS1 1 


0 

0 

0 

0 

0 

0 

23 

S2*IS2l 


0 

-8.325298e-07 

0 

+ 1.574352e-06 

0 

0 

24 

RM*lRMl 

a i 2 = 

-8.609022e-08 

0 

— 1 .771 749e— 07 

-1.772552e-07 

— 4.532282e-08 

0 

25 

AF*lAFl 

0 

0 

0 

0 

0 

0 

26 

N1*N2 


0 

-1.739066e-06 

0 

0 

0 

0 

27 

N1 *S1 


0 

0 

0 

0 

0 

0 

28 

N1 *S2 


0 

0 

0 

0 

0 

0 

29 

N1*RM 


0 

0 

0 

0 

0 

0 

30 

N1 *AF 


0 

0 

0 

0 

0 

0 

31 

N2*S1 


0 

0 

0 

0 

0 

0 

32 

N2*S2 


0 

0 

0 

0 

0 

0 

33 

N2*RM 


0 

0 

0 

0 

0 

0 

34 

N2*AF 


0 

0 

0 

0 

0 

0 

35 

S1*S2 


0 

0 

0 

0 

0 

0 

36 

SI *RM 


0 

0 

0 

0 

0 

0 

37 

SI *AF 


0 

0 

0 

0 

0 

0 

38 

S2*RM 


0 

0 

0 

0 

0 

0 

39 

S2*AF 


0 

0 

0 

0 

0 

0 

40 

RM*AF 


0 

0 

0 

0 

0 

0 

41 

INUN2I 


0 

0 

0 

0 

0 

0 

42 

IN1+S1 1 


0 

0 

0 

0 

0 

0 

43 

1 N 1 *S2 1 


0 

0 

0 

0 

0 

0 

44 

iNURMl 


0 

0 

0 

0 

0 

0 

45 

IN1+AFI 


0 

0 

0 

0 

0 

0 

46 

lN2*S1 1 


0 

0 

0 

0 

0 

0 

47 

lN2*S2l 


0 

0 

0 

0 

0 

0 

48 

lN2*RMl 


0 

0 

0 

0 

0 

0 

49 

lN2*AFl 


0 

0 

0 

0 

0 

0 

50 

IS 1 *S2l 


0 

0 

0 

0 

0 

0 



(C 

>nly 50 of 

37 regression coefficient rows shown) 



Fig. 3b Iterative Method: Coefficients of optimized regression models of gage outputs Rl, R2, R6. 
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MATH 

TERM 

T— STATISTIC 
(from Fig. 2a) 

EXACT SOLUTION 

(uses coefficients of non-iterative method, 
see also Eq. (4), values given in Fig. 2b) 

APPROXIMATION 
(uses coefficients of iterative method, 
see also Eq. (9), values given in Fig. 3b) 

STD. DEV. OF 
LOAD RESID., % 

INTERCEPT 

+244 

770 = +5.650# + 01 

~ a ° = +5.663#+ 01 
a 1 

0.0489 % 

Rl 

+982 

771 = +2.915# + 00 

— = +2.909# + 00 
a 1 

0.0566 % 

R2 

+29 

772 = +7.997#- 02 

~ a , 2 = +8.065# 02 

a 1 0 2 

0.0491 % 

i?3 

+4 

773 = +2.409# -03 

~ a 3 = +3.105# 03 

a 1 c 3 

0.0492 % 

#5 

-11 

774 , = -1.742# -02 

~ U4 = 1.838# 02 

a 1 e4 

0.0495 % 

m 

-27 

775 = -2.561# -02 

. = 2 . 444 # 02 

a 1 |«i| 

0.0492 % 

\R2\ 

-18 

77s = -1.567# -02 

~ a 7 6 . = 1.395# 02 

a 1 \b 2 \ 

0.0496 % 

|i?.3| 

+13 

777 = +8.287# -03 

= +7.948# 03 

ai |c 3 1 

0.0490 % 

\R5\ 

+14 

77s = +2.498# -02 

_a8 , = +2.336# 02 

a 1 e 4 

0.0502 % 

R5-R5 

-5 

t? 9 = -8.925# -06 

~ ag = 8.157# 06 

cl 1 e4 e4 

0.0492 % 

Rl ■ |#1| 

-7 

7710 = -3.293# -05 

_ai ° = 3.747# 05 

ai ai \a\ \ 

0.0506 % 

R2 ■ \R2\ 

-8 

7711 = -3.187# -05 

~ ai \ . = 3.257# 05 

ai 62 p 2 | 

0.0490 % 

R5 ■ |i?5| 

+5 

7712 = +8.664# — 06 

_ai , 2 . - +9.327# 06 

ai e 4 e 4 

0.0492 % 


Fig. 4 Comparison of exact solution with approximation of regression model coefficients of N 1. 
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